A toy model for molecular condensates in Bose gases 
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The occurrence of a molecular Bose-Einstein condensate 
is studied for an atomic system near a zero energy resonance 
of the binary scattering process, with a large and positive 
scattering length. The interaction potential is modeled by a 
pseudo-potential having one bound state. Using a variational 
Gaussian ansatz for the iV-body density operator, we discuss 
the thermodynamic properties at low temperature and the 
relative stability of the system towards the formation of an 
atomic Bose-Einstein condensate. We also derive an approxi- 
mate Gross-Pitaevskii equation for the molecular condensate 
leading to the prediction of a Bogoliubov spectrum. 

The atomic Bose-Einstein condensation has been ex- 
tensively studied during past years both theoretically and 
experimentally [0J^]. In this Letter, we consider the dif- 
ferent case of a molecular Bose-Einstein condensate. The 
realization of such a situation is a challenge in the field of 
low temperature physics. Already cold gases of molecules 
have been produced |^-[|. No evidence for a molecular 
condensate has been obtained yet, although the stimu- 
lated Raman technique used in j^], forming molecules 
directly from an atomic condensate, is promising. 

Another scenario is possible for the formation of the 
molecular condensate. Consider the case of a so-called 
zero-energy resonance in the two-body scattering process, 
where the positive scattering length a > is much larger 
than the effective range r e of the interaction potential 
In this case the two-body potential supports a s- 
wave bound state </>o with a spatial extension ~ a much 
larger than the other s-wave bound states and with a 
much smaller binding energy Eq — — fi 2 /(ma 2 ). One then 
would rely simply on three-body collisions between atoms 
of a trapped atomic condensate to produce molecules in 
the state 4>q. As shown in the rate of formation of 
molecules scales as a in the limit of large a and mainly 
leads to formation of molecules in the highest s-wave two- 
body channel Q that is in the state 4>q. From energy 
conservation one finds that the molecules produced in 
4>q have a kinetic energy —Eq/3 when the initial atomic 
wavevectors k are such that fca < 1. The requirement 
that the molecules formed in state cj>o remain in the trap 
imposes a depth of the trapping potential larger than 
—Eq. For a modest trap depth of 10 fiK achievable in an 
optical trap || and for the mass of Rubidium one finds 
that a has to be larger than 445 Bohr radii. Such high 



values of a may be obtained using a Feshbach resonance 
induced by a magnetic field |l(||ljfl . 

In this Letter we assume that a condensate of molecules 
in the diatomic bound state 4>q has been formed. We de- 
scribe the molecular condensate at the atomic level: we 
take as a starting point a model Hamiltonian for interact- 
ing atoms and we use a Gaussian variational ansatz for 
the ./V— body density operator including only the effect 
of binary interactions between atoms. In the low den- 
sity regime (na 3 <C 1 where n is the density of atoms) 
we determine (i) the properties of the molecular conden- 
sate at thermal equilibrium and (ii) the response of the 
molecular condensate to a small perturbation of the trap- 
ping potential. Our results for both cases correspond 
formally to the known properties of a weakly interacting 
Bose-Einstein condensate of particles of mass 2m having 
a coupling constant 6g, where m is the atomic mass and 
g = AirTi 2 a/m is the atomic coupling constant. At higher 
density (na 3 > 7r/192) we find that an atomic condensate 
forms. To model the binary atomic interaction potential 
we use the pseudo-potential V defined by the following 
action on a two-body atomic wavefunction i/j: 

(ri,r 2 \V\ip) =gS(f)d r \rip(R - f/2, R + r/2)] , (1) 

where we have introduced the coordinates of the center 
of mass and of the relative motion: 



5 n+r 2 _ _ _ 

R = , r = t*i — ri and r = r 
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This model potential discussed in jl3] has been used to 
extend the BCS theory to inhomogeneous atomic sys- 
tems As we now see, it has the ability to capture 
the essential feature of a general binary scattering prob- 
lem near a zero-energy resonance. Consider indeed the 
relative motion of two atoms described by the hamilto- 
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Contrary to the case of the usual contact interaction, the 
scattering of two atoms interacting with the potential 
(jl|) is a well defined problem. The scattering amplitude 
for a relative wavevector k is given by — a /{I + ika) re- 
producing the known universal Lorentzian shape of the 
scattering cross section near a zero energy resonance [pj . 
In addition to the usual scattering states, the pseudo- 
potential for a > leads to the asymptotic form of the 
s-wave bound state of energy Eq: 
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This property of the pseudo-potential reproduces the uni- 
versal behaviour for a general interaction potential close 
to a zero-energy resonance, already mentioned in the in- 
troduction. It plays an essential role in our approach as 
it allows the formation of pairs of atoms, that is binary 
molecules, while keeping the mathematical simplicity of 
a zero range model potential (r e = 0). 

The goal of this paper is to study the properties of a 
molecular condensate, in a symmetry breaking approach, 
this corresponds to = and (4^) ^ 0, where ^> is the 
atomic field operator. However it is dangerous to exclude 
a priori the coexistence of an atomic and a molecular 
condensate. We therefore use the more general symmetry 
breaking prescription by splitting the field operator in a 
classical part and a quantum fluctuation part: ^ = $+</>, 
with (4>) = 0. We choose a Gaussian ansatz for the many- 
body density matrix D = Z^ 1 exp —j3K with 
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+ - I d 3 r 1 d 3 r 2 [4>(r 1 )A(r u r 2 )ct>(r 2 ) + h.c] . 

There are three variational fields in this theory: h, A 
and $. At thermal equilibrium, the grand potential 
E — jiN — TS is minimum so that the functions h(fi, f 2 ) 
and A(r*i, f 2 ) may be expressed in terms of the one-body 
density matrix p(fi, r 2 ) = (r 2 )(f)(ri)) and of the pair- 
ing function R,(fi,f 2 ) = (4>(ri)(j)(f 2 )}, while <J>(r) verifies 
a partial differential equation p4j 

h{r u r 2 ) = _^l(V 2 5)(r ) + (2gn(R) - p)S(r) (5) 
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In these equations, we have introduced the atomic density 
at a point R 



n(R) = \<£\ 2 (R) + p(R,R) , 
and the regular part of the pairing function 

R IC JR ) = lim d r \rR{R - r/2, R + r/2) 

r—>0 L 



(8) 
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We assume here that the atoms are in a cubic box of 
size L with periodic boundary conditions ifTBf , so that 
$, n, K reg do not depend on position. We expand the field 
operator on plane waves using a Bogoliubov transform: 
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The commutation relations of the bosonic annihilation 
operators 6^'s lead to the normalization of the modes 
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1. We search the (uk,Vk)''$ so 



that if is a sum of decoupled harmonic oscillators: 

K = Ko + J2 nuJ kbtb j ; . (11) 

k 

From the equilibrium conditions (|^,^), we find that each 
(itfc, is the eigenvector of the system: 
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with the notation hi. 
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2gn — //. A simple al- 
gebra gives the expressions for the eigenvectors and the 
spectrum {htok} |l6[: 
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huj k = [h 2 k ~ g 2 ($ 
The pairing function is given by 



(13) 
(14) 
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and the one-body density matrix 



p(n, r 2 ) = ^ E + X ) v " + M ex P(^ • ■ (16) 



In Eqs.( |15yiq ), fk is the Bose occupation factor fk = 

[exp(/37iw fc ) - l]" 1 . 

All the equilibrium properties may be expressed in 
terms of (n,T). For this purpose, we have to determine 
the three unknown parameters (/i, $,K reg ). In the pres- 
ence of an atomic condensate, $ does not vanish and a 
first relation is given by Eq.(^) 



/i = g[2n + K r 



A second one is obtained from Eq. 
in Eq.© 
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by setting r = 
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The third equation is obtained by extracting the regular 
part of Eq. ([L5j) as in |13| ; this leads to 
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In what follows, we consider the thermodynamical limit 
so that the sums over k are replaced by integrals, we also 
restrict to the case of a vanishing temperature T = 0, so 
that the Bose occupation factors fk is zero. 

First, we note that the stability of the ground state 
imposes real values for the spectrum {Tkjfc}, this implies 
K rcg < and from Eq.(|l9|), we find |re re g| > \^\ 2 - In par- 
ticular for a vanishing number of atoms in the condensate 
$ = 0, we get from Eqs.(0g|||): 

(20) 
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This value of the density determines the threshold of co- 
existence of an atomic condensate (<£> ^ 0) with a molec- 
ular one. For n > n c , the model predicts the coexistence, 
a result already obtained in a slightly different approach 
in [p"7f . In this high density regime, our modelization 
is questionable and we do not pursue the study in this 
range anymore. 

For n < n c , there is no atomic condensate: $ is iden- 
tically zero so that Eq.(fl7|) does not hold in this regime 
and equilibrium properties are deduced from Eqs.([l8 19) 
only. We now wish to check that the low density regime 
corresponds indeed to a molecular condensate. For this 
purpose, we suppose that the chemical potential tends 
to a finite negative value, hence 3> gn,g\R re g\ and 
huj k ~ Ti 2 k 2 /{2m) - fi. From Eq.(|l|) we find the lowest 
order approximation 
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which is finite indeed, and negative. This result is en- 
lightening: Eq is just the binding energy of a pair of 
atoms in the bound state Eq.(Q) so that the value Eq.(|2l|) 
of the chemical potential corresponds to that of an ideal 
gas of molecules. Since we are at zero temperature this 
gas of pairs is actually a molecular Bose-Einstein conden- 
sate. 

This interpretation is confirmed by the lowest order ap- 
proximation to the pairing function n in the low density 
limit. One first calculates R rcg from Eq.(|l8|) using the 
lowest order approximation Vk — gS reg (-Eo — h 2 k 2 /m) _1 
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then we calculate the integral over k in Eq.(|l5|) and in 
Eq.(|l6|) to obtain the lowest order contributions to the 
pairing function and to the one-body density matrix: 



K(ri,fa) ^ \fn <f> (r) , p(fi, f 2 ) ~ nexp (--J 



(23) 



This expression of k(t*l,7%) clearly shows that the pair- 
ing function describe the spatial structure of two atoms 



linked in the molecular bound state <pQ. Similar results 
have been obtained in jl7j for a different model potential. 

For a small but finite gaseous parameter na 3 the molec- 
ular condensate is not an ideal gas but rather a weakly 
interacting Bose gas, with an effective coupling constant 
ffmoi between the molecules. We derive this coupling con- 
stant by calculating the first correction to the expression 
Eq.(pTl) for the chemical potential. We expand the ex- 
pression Eq.(H4) to first order in na 3 and substitute the 
result in Eq. (jl9|) . This leads to 



fi= -y(l - 127ma 3 + . . .) 



(24) 



On the other hand the molecular chemical potential ^ mo i, 
equal to twice the atomic chemical potential is given 
in the usual mean-field approach for condensates by 



Mmol — ^0 + 5mol n mol 



(25) 



where n mo i = n/2 is the density of molecules. We there- 
fore deduce for the coupling constant between molecules 



5moi = 6g 



(26) 



As (? mo i is positive the molecular condensate is stable with 
respect to a spatial collapse. 

The value of the molecular coupling constant g mo i can 
also be obtained from the response of the gas to a time 
dependent perturbation. It is important to check that 
the corresponding value coincides with the static pre- 
diction Eq.(p6|). Imagine that one perturbs the system 
from thermodynamical equilibrium by applying an exter- 
nal potential on the atoms for a finite time interval [0, r]. 
We describe the evolution of the gas by a time dependent 
Gaussian ansatz for the many-body density operator [fl4"| . 
In that case, as in jljj, we deduce from the Heisenberg 
equation for the field operator and from Wick's theorem 
the time evolution of the pairing function, written here 
for convenience for t > r: 
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A n + 2(gn(n) + gn{f 2 ) -fj,)+H, 



?K reg (fi)p*(f 1 ,f 2 ) -3K rcg (r 2 )p(r 1 ,r 2 ) . (27) 



We assume that the applied perturbation varies very 
slowly spatially at the scale of the scattering length a. 
Thus the pairing function has negligible components on 
the scattering states and can be assumed to have the 
same r dependence as 4>o- We therefore take the Local 
Density Approximation: 



1/2 

K(ri,r2)= n{R,t) 4>o{r) exp iS(R,t) 



(28) 



We close equation (g7|) with the Local Density Approxi- 
mation for the one-body density matrix: 
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p(ri,r 2 ) = n(R,t) exp 



(29) 



A simple projection of Eq.(p7|) on the bound state (Q) 
leads to the Gross-Pitaevskii equation 



ihdtipp = 
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A A + 6g\il> P \ J -2n + Eo)il)P , (30) 



where we have introduced the macroscopic wave func- 
tion ipp{R,t) = (n/2) 1 exp(zS') describing the molec- 
ular condensate. We note that this equation, confirms 
in a direct way our previous finding on /j, at equilibrium 
(Eq.(pij)). From the linear analysis of Eq.(|30"|), we pre- 
dict a Bogoliubov spectrum for the molecular condensate 
different from the atomic condensate: 
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Measurement of this spectrum could be used as an ex- 
perimental evidence for the condensation of pairs. 

As a conclusion let us stress three points. First, one 
word about the temperature. Indeed, the critical tem- 
perature at which the condensate of molecules (mass 2m, 
density n/2) forms is 



k B T c = 

m 



2C(3/2) 
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(32) 



In the low density regime kpT c <C |£b|, so that condensa- 
tion of pairs can occur without any thermal dissociation. 
Second, it would be interesting to test the prediction on 
the coupling constant (6g) by a direct calculation of the 
scattering of two molecules of the type considered here. 
Finally, our model does not describe the relative stabil- 
ity of this molecular condensate toward the formation 
of molecules in deeper bound levels. An evaluation of 
the creation rate of deep bound states by collision of an 
atom with one molecule would be a relevant complement 
to this analysis. 
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tory. 
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